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Abstract 

We report the existence, and study mobility and interactions of gap polariton solitons in a 
microcavity with a periodic potential, where the light field is strongly coupled to excitons. Gap 
solitons are formed due to the interplay between the repulsive exciton-exciton interaction and 
cavity dispersion. The analysis is carried out in an analytical form, using the coupled-mode (CM) 
approximation, and also by means of numerical methods. 
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I. INTRODUCTION 



The strong light-matter coupling in semiconductor microcavities has recently attracted 
much attention [Tj. In particular, the strong and fast nonlinear response of microcavity 
exciton-polaritons has allowed to predict and observe several important nonlinear effects, 
such as bistability [21 El H], parametric wave mixing [3j IH [5], superfluidity (H E] and for- 
mation of solitons [8j [9]. While the exciton-polariton nonlinearity is defocusing due to 
the electrostatic repulsion of excitons, the effective dispersion of the electromagnetic wave 
may be controlled in microcavities with periodic potentials. The latter can be created by 
mirror patterning [TD], or by way of surface acoustic waves [TTl IT2"] . In either case, the peri- 
odic modulation of system parameters can be achieved on the micron scale, leading to the 
emergence of gaps in the polariton spectrum. Localized nonlinear modes with the Fourier 
transform residing within the forbidden gaps of linear spectra are called gap solitons (GSs), 
also known as Bragg solitons. GSs may exist with any sign of the nonlinearity. Photonic GSs 
have been extensively studied in fiber gratings and planar optical lattices [13] . Matter- wave 
GSs have been observed in the atomic condensate of 87 Rb [13]. From the vast literature on 
solitons in periodic structures it is relevant here to mention works which either considered 
cavity effects or where material excitations played a crucial role. These include the soliton 
transmission through resonantly absorbing Bragg gratings [T5j \1Q\ UT] , control of electro- 
magnetically induced transparency using photonic bandgaps [18J, and light-only solitons in 
microcavities with photonic crystals [T9| [2D] . The aim of this work is to initiate studies 
of the exciton-polariton GSs, which are half-light half-matter nonlinear excitations, whose 
self- localization is supported by a periodic potential acting on the photonic component. 



II. THE POLARITON MODEL AND ITS LINEAR PROPERTIES 

Below we focus on the microcavity model, which disregards cavity losses, aiming at the 
proof-of-the-principle demonstration of the existence and robustness of gap polariton solitons 
in this setting. Effects of the dissipation and introduction of a compensating gain may be 
important to the experimental realization, and will be considered elsewhere. In the scaled 
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form, the equations for local amplitudes of the photon (E) and exciton (ty) fields are [TJ [7] 

8 t E -i{8l + d 2 y )E - iU(x, y)E = M , (1) 
d t ^ +i\m\ 2 m = %E . (2) 

In these equations, the time and coordinates x, y are measured, respectively, in units of 
1/VLr and j '(2Qr), where £Ir is the Rabi frequency, u and k = cn/oo are the pump 
frequency and wavenumber, with n being the refractive index. Further, Qn\E\ 2 /g and 
Qn\^\ 2 /g are numbers of photons and excitons per unit area [5], and g is the exciton- 
exciton interaction constant. Taking typical parameters of a microcavity based on a single 
InGaAs/GaAs quantum well, hVtR ~ 2.5 meV, hg ~ 10~ 4 eV-/im 2 (2J Ej, one finds that 
\E\ 2 = 1 corresponds to the electromagnetic field with intensity ~ 10 kW/cm 2 , while the 
time and length units translate into ~ 0.6 ps and ~ 1 /im, respectively [8], |9]. 

The form of Eqs. Q and Q, with zero detuning between them, assumes identical reso- 
nance frequencies of photons and excitons. Note that the separation between the nonlinearity 
and diffraction in these equations resembles the phenomenological model introduced earlier 
in Ref. [21]. 

The last term on the left-hand side of Eq. Q is the lattice potential induced by the 
periodic modulation of the cavity resonance [10]. First, we consider one-dimensional (ID) 
case, with ^/-independent fields, and take potential 

U(x) = ecos(2/c x), (3) 

where e is the depth of the potential and 7r//co is its period, so that the first Brillouin zone 
for polariton momentum k is — k < k < k . A 2D model, in which U(x,y) is periodic in x 
and localized in the y direction, is considered towards the end of the paper. 

Without the lattice potential, e = 0, solutions to the linearized version of Eqs. and 
2|) are sought for as E, ^ ~ ^kx-iAt^ yields the spectrum consisting of two branches, 



A ± (*) = ^±^ + l, (4) 

see Fig. [TJa). The addition of the lattice potential splits this spectrum into multiple bands 
with the zone folding happening at k — fco, see Fig. []Jb). Gaps between the bands are 
getting wider for deeper lattices (larger e). Unusually, the choice of k Q also affects the gap 
widths, see Fig. [l](c) . This happens because the curvature of the exciton-polariton dispersion 



FIG. 1: (Color online) (a) The polariton spectrum of the homogeneous microcavity (e = 0). Black 
(red) line corresponds to the lower (upper) polariton branch, (b) The spectrum of the microcavity 
equipped with the periodic potential: e = 1, ko = 7r/3. (c) fco-dependence of the width of the 
principal gaps in the lower (black solid lines) and upper (dashed red lines) polariton branches, for 
e = 1. (d) The polariton spectrum in a vicinity of the primary gap of the lower polariton branch. 
Solid black lines show the the exact spectrum, and dashed red lines represent the approximation 
obtained from Eqs. (|6]), 0. 

without the periodic potential, see Eq. Q, strongly depends on k. In contrast, the photonic 
dispersion is parabolic, and it is modified by the periodic potential in such a way that the 
width of the primary gap is A^ = e, being obviously independent of k . In Fig. [TVc) we plot 
the widths of the primary gaps in the two dispersion branches (A±) as functions of ko- For 
relatively small ko, the gaps in the upper and lower polariton branches are approximately 
the same. Increasing ko, the dispersion of the upper polariton branch tends to its photonic 
(light-only) limit, hence the width of the primary gap in the upper branch increases and 
tends to A^ = e. For the lower polariton branch, which is the practically important one 
(see below), A^ first increases and then drops to zero for large ko, see Fig. []Jc). Below we 
focus on GSs residing in the principal gap on the lower polariton branch. Our choice of ko 
throughout this paper is tt/3, which corresponds to modulation period ~ 3/xm and matches 
the experimental conditions of Ref. |10j . 

III. THE EXISTENCE AND ROBUSTNESS OF GAP SOLITONS IN THE ID 
GEOMETRY 

In order to address the existence of GSs in the model, we first apply the coupled-mode 
(CM) approximation, which is known to lead to explicit analytical results [13J, which are 
then compared to full numerical solutions of Eqs. (IT]) , ((2]) . To this end, we introduce 



F = [E, \I/] T and assume 



F = cio [C + (x, t) exp(— iA$t + ikox)+ 
C_(x, t) exp(— iA t — ik x)} , 



(5) 



where A = A_(A; ) is the frequency as given by Eq. Q in the middle of the gap, 
do = [— Ao, 1] T is the polariton eigenvector in the absence of the nonlinearity and peri- 
odic potential, and C + and C_ are slowly varying functions. The substitution of Eq. ^ 
into Eqs. Q, ^ and manipulations similar to those known in the context of the CM equa- 
tions for fiber Bragg gratings [131 123], we derive the CM equations corresponding to the 
present setting: 



i (A 2 , + 1) d t C+ + 2ik A 2 d x C + + kC„ 



|C 4 



2|C_| 2 ) C 4 



i (Aq + 1) d t C- - 2ik A 2 d x C- + kC^ 
- (|C_| 2 + 2 |C+| 2 ) C_ = 0, 



(6) 



(7) 



where k = cAq/2 is the coefficient of the Bragg-reflection-induced linear coupling between 
the counterpropagating waves. Figure [Ud) compares the spectra found from the full model 
and from the CM equations. The good agreement between the two persists for ejk\ < 1, 
i.e., for relatively weak potentials. 

Using variables £ = x/ (2/coAq) and t = t/(l + Aq), we obtain from Eqs. ^ and (J7| 
explicit solutions for GSs [T3] . 

C± = u±{rf) exp iqij]) ± is(r]) — i5\/l — V 2 r , (8) 



u 



2[k 2 -5 2 )(l + V)VT^V 2 ~ 
(3 - V 2 ) /tcosh(2 v / ^ 2 - P?7) - S 



l-V 2 



— tan 



q = V5r] + 



AV 



3-V 2 



K + 5 
K — 5 



tanh {^J k 2 — 5 2 r] 



(9) 
(10) 

(11) 



where 5 is the frequency detuning relative to the gap center, \V\ < 1 is the soliton velocity, 



and V = (^~ Vt)/VI~V 2 . 
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FIG. 2: (Color online) Numerical (solid black lines) and analytical (dashed red lines) solutions for 
gap solitons at e = 1, ko = tt/3, A = —0.7. 



In Fig. |2| we compare the analytical stationary solitons given by Eqs. (|8j) and their 
counterparts found numerically from Eqs. ([I]) and (p2J). The overall agreement is reasonable. 
The main source of the error is the inaccuracy of ratio which, in the framework 

of the CM approximation, is taken as per the linear eigenvector at e = and A = Ao, 
therefore it is assumed to remain constant while the soliton's spectrum is shifting within the 
gap, following a variation of 5. However, numerical solutions show a tangible dependence of 
the l^l/l^l ratio on the soliton frequency 5. 

Next, we check the stability of the numerically found GS solutions. To this end, we per- 

rp rp 

turb them by setting [E, = [E Q (x) + e(x,t), ^q(x) + ^(x,t)\ -exp(— iA t), and linearize 
Eqs. ([I]), ([2]) assuming that perturbations e (x, t) and ip { x ,t) are small. This yields 

(12) 
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where L e = A + d% + ecos(2A;ox). Then, we seek solutions to Eq. (12) as y(x,t) = 
exp(\t)y (x) . The existence of eigenvalues with Re {A} > implies instability. 

We have found that results of the stability analysis for the full polariton model qualita- 
tively coincide with the known stability properties of the GSs obeying CM equations (|6]), 
([7]) [221 In particular, the gap polariton solitons are stable in the lower half of the gap, 
and feature various instabilities in the upper half, see Fig. [3|a). Figure [3](b) demonstrates 
that the instability (if any) initiated by random perturbations causes the soliton to ramble 
erratically across the lattice. 

We have also checked numerically mobility and collisions of the GSs in the full model. 




Re(K) t 

FIG. 3: (a) The stability spectrum for the soliton with A = —0.5 and other parameters as in Fig. [2} 
Eigenvalues with Re {A} > correspond to the instability, (b) Evolution of an unstable soliton 
perturbed by random noise. 
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FIG. 4: (a) Pinning of the soliton, which was initially moving with a relatively small velocity, 
V = 0.1 (b) Merger of colliding in-phase solitons moving with velocities V = ±0.2. (c) Rebound 
of 7r-out-of-phase solitons colliding with the same velocities as in (b). (d) Collision between the 
soliton moving at V = 0.2 and a quiescent one. Parameters of the potential are the same as in 
Fig. 1(b), (d), and A = -0.6. 

Figure |4^a) shows a soliton which initially moves through the lattice as prescribed by the 
CM approximation, but eventually gets pinned around one of the lattice sites. Outcomes of 
collisions between the solitons are sensitive to both initial velocities and the relative phase of 
the interacting solitons, as shown in Figs. |4](b)-(d) . Collisions between solitons with opposite 
velocities result in merger of in-phase soliton pairs, and rebound of the solitons with the 
phase different of tt, see Figs. |4]^b)-(c). Collisions of solitons with different velocities can 
produce a plethora of outcomes, with one example shown in Fig. |4|c). 

IV. GAP SOLITONS IN THE 2D GEOMETRY 

To examine the relevance of the ID model elaborated above to the 2D geometry of 
practical interest, we consider a configuration where the periodic potential acting in the x 
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FIG. 5: (a) The surface plot of the potential (13) used in the 2D model, for e = 1, ko = 0.85, 



w = 2 (the usual notation of the quantum theory requires to replace U in Eq. ([I]) by —U, so that 
the potential structure in (a) corresponds to the periodic sequence of potential wells), (b) The 
soliton's profile for A = —0.65. 

direction is combined with a localized potential applied along the y coordinate: 

U{x, y) = -e{l- exp [-{y/wf] [1 + cos(2A;ox)] } . (13) 

The shape of the above potential is shown in Fig. [5^ a). We have found the corresponding 
2D soliton solutions numerically, see Fig. |5](b), using a time-independent iteration method. 

To check if the dynamics seen in the ID case is retained in the 2D configuration we have 
carried out a series of numerical experiments in soliton collisions. The initial conditions were 
set as 

E = = Ae- {x - xo)2/w *- y2/w2 y cos{k x)e ikx , (14) 

where A, w x and w y have been chosen to approximate the stationary profile of the numeri- 
cally found solitons and k is the initial soliton momentum. The results of these simulations 
are shown in Figs. [6] and [7} Similar to the ID case, the in-phase solitons tend to merge in 
the course of the collision, while the out-of-phase solitons bounce back. The outcome of the 
collision of the out-of-phase 2D solitons is similar to what was observed in the ID model. 
However, for the in-phase solitons the large-amplitude pattern generated by the merger is, 
most often, unstable, splitting into two quasi-solitons moving in opposite directions, see 
Fig. [6j This pair is asymmetric, one of the emerging quasi-solitons being usually larger and 
slower than the other. 



V. SUMMARY 

We have predicted the existence and studied stability, mobility and interactions of gap 
polariton solitons in microcavities equipped with periodic photonic potentials. The ID 
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FIG. 6: Collision of in-phase solitons in the 2D geometry. All panels show \E(x,y,t\ 2 . The top 
panel shows dynamics in the (x,i)-plane along y = 0. Bottom panels show the light intensity in 
the (x, y)-plane for fixed time. Parameters are the same as in Fig. [5] Initial conditions for each 
soliton are chosen as in Eq. (14), with A = 0.5, w x = 5, w y = 3, k = ±ko/20. 




FIG. 7: The same as in Fig. |6j but for the it out-of-phase solitons. 



model has been studied using the CM (coupled-mode) approach, which yields analytical 
solutions for the solitons, and also by way of the numerical solution of the full system of 
equations for the photonic and excitonic components of polaritons. Furthermore, we have 
studied the two-dimensional microcavity with the periodic potential along one dimension 
and the trapping potential along the other. 
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